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ABSTRACT 

We report integral-field spectroscopic observations with the Cambridge Infra-Red Panoramic 
Survey Spectrograph (CIRPASS) mounted on the GEMINI South telescope of the nucleus of 
the nearby galaxy NGC 5128 (Centaurus A). We detect two-dimensional distributions of the 
following emission-lines: [PII], [Fell] and Paschen (3. We compare our observations with pre- 
viously published radio observations (VLA) and archival space-based near-infrared imaging 
(HST/NICMOS) and find similar features, as well as a region of high continuum coincid- 
ing with the jet (and its Nl knot) at about 2" North-East of the nucleus, possibly related to 
jet-induced star formation. We use the [FeII]/[PII] ratio to probe the ionisation mechanism, 
which suggests that with increasing radius shocks play an increasingly important role. We 
extract spatially resolved 2D kinematics of Pa/3 and [Fell] emission-lines. All emission-line 
regions are part of the same kinematic structure which shows a twist in the zero-velocity curve 
beyond ~ 1" (for both Pa/3 and [Fell]). The kinematics of the two emission-lines are similar, 
but the Pa/3 velocity gradient is steeper in the centre while the velocity dispersion is low ev- 
erywhere. The velocity dispersion of the [Fell] emission is relatively high featuring a plateau, 
approximately oriented in the same way as the central part of the warped disk. We use 2D 
kinematic information to test the hypothesis that the ionised gas is distributed in a circularly 
rotating disk. Assuming simple disk geometry we estimate the mass of the central black hole 
using Pa/3 kinematics, which is consistent with being distributed in a circularly rotating disk. 
We obtain M. = 8.25±|;|| x 10 7 M , for PA = -3° and i = 25°, excluding the M. - a 
relation prediction at a 3<r confidence level, which is in good agreement with previous studies. 

Key words: galaxies: elliptical and lenticular - galaxies: kinematics and dynamics - galaxies: 
structure, galaxies: individual, NGC 5128 



1 INTRODUCTION 

NGC 5128 is the king galaxy of the southern hemisphere, if 
sovereignty is judged by the number of studies devoted. It is no 
wonder that this is the case since Centaurus A, or Cen A, as it 
is co mmonly referred to, is the fifth brightest galaxy in the sky 
(RC2, de Vauc ouleurs et al. usually classi fied as an EOp or 

SOp. It hosts a powerfu l X-ray and radio jets ( Hardcast le et al 



20031: lEvans et al.ll2004l) and a low-luminosity AGN jjones et al 
1996l :lchiaberge et alj200l"h . Cen A suffered a recent major merger 



that shaped its pr esent appearance, leaving shells at large radii 
(Malin et al. 198 3]) and a conspicuous polar dust lane in the cen- 



tre dDufour et alj[l979l). consistent with being a thin warped disk 
dOuillen et al.l 1 19921 ; iNicholson et all Il992l : louillen et alj [19931 
2006). It is, however, the proximity that makes this object truly spe- 
cial. Its most recent d istance estimate places it at 3.42 ± 0.25 Mpc 
dFerrarese et al.l2 006). while other recent studies give : 3.40 ± 0.15 
Mpc dlsraelll 19981). 4.2 ± 0.3M pc dTonrv et alJboOll) and 3.84 ± 
0.35 Mpc dReikuba et alj|2005h . Being so close, Cen A is also the 
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closest AGN and recent merger, a great case for testing our under- 
standing of processes that shape galaxy formation and ev olution 
(for d etailed information about the galaxy see the review bv llsraell 
Il998h . 

The thick dust lane that crosses Cen A, making it a show case 
for astronomical PR, hides the centre of the galaxy an d its con- 
stituents. Within the AGN paradigm dAntonuccil 1993), it is ex- 
pected that there is a suppermassive black hole at the bottom of 
Cen A's potential well. Indeed, recent studies, 'peering through 
the dust' with near-infrared spectrographs, investigated this as- 
sumption and found that it is necessary to place a dark object in 
the nucleus with a mass of: M. = (2±f°) x 1 8 Mp) for in- 
clination i > 15° (ground-based gas kinematics, iMarconi et al.l 
|2001|) . M. = (2.4±g;|) x 10 8 M Q for an edge-on model, M. = 
(1.8±o.4) x 10 8 M Q for i = 45° and M. = ( LjSlgj) x 10 8 M B 
for i = 20° (ground-based stellar kinematics, Silge et alj|2005l) . 
M. = (6.lto j) x 10 7 M Q at i = 45° (gr ound-based adaptive op- 
tic as sisted observations of gas kinematics, Harin g-Neumaver et al.l 
120061) and M. = (1.1 ± 0.1) x 1O 8 M for i = 25° and 
M. = (5.5t° Q 7 6 ) x 1O 7 M for i = 45° (ground- and s pace- 
based observations of gas kinematics. IMarconi et al.ll2006j) . The 
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last study also constrained the size of the central massive object 
to r ~ 0.6 pc suggesting that it is indeed a black hole. The im- 
portance of a robust measurement of Cen A's M. is highlighted in 
the fact that the above masses place the black hole in Cen A ap- 
proximately betw een 2 times and an order of magnitude above the 
M. - a relation dGebhardt et al.ll2000l: IFerrarese & Memtdl200(j| ; 



Table 1. Observational details 



iTremaine et~al]|2002l ; IFerrarese & FordH2005h . Although the large 
scatter in the relation and the low number statistics might be mis- 
leading, i t is also possible that Cen A, as a recent merger, is a spe- 
cial case dSilge et al.ll2005h . 

Each of these studies brought an improvement in spatial and 
spectral resolution, but they were all based on observations with 
long-slits along a few position angles. For construction of dynam- 
ical models, whether stellar or gaseous, it is crucial to have suffi- 
cient information about the behaviour of the constraints. In the case 
of stellar dynamical models, it is necessary to observe the kinemat- 
ics at all pos itions (within ~ lR e ) to robustly constrain the orbital 
distribution dCappellari & McDermidll2005l ; iKrainovic et alj|2005t 
Shapiro et al. 200ej), while for gaseous models the geometry of the 
gas disk can only be confirmed by two-dimensional mapping and 
only a well-orde red gas disk can be used to reliably determine the 
enclosed mass flSarzi et al.l l2002). 

The nucleus of Cen A is active and the dusty and gaseous po- 
lar disk is significantly warped from small to large scales, where the 
disk alternates between having the southern side (within 60") and 
northern side (b e yond 100") nearest t o the observer (see F ig. 1 1 of 
iNicholson etai] Jl992h and Fig. 5 of lOuillen etall d2006l) ). When 
constructing dynamical models it is important to understand in 
what way these two properties shape the observations. The ioni- 
sation of emission-lines could be under a strong influence of the 
AGN, while the kinematic signature of a disk warp could be hard 
to recognise from long-slit measurements only. 

In this paper we present the first observations with an integral- 
field spectrograph of the nucleus of Cen A in the near-infrared. In 
Section [2] we describe the observations and the data reduction. In 
Section[3]we analyse the data and compare them with existing ob- 
servations. We describe the kinematical properties of the nuclear 
disk-like structure in Section [4] In Section [5] we discuss the evi- 
dences for the central black hole. Section [6] summarises our con- 
clusions. 

Throughout the study we assumed a distance to Cen A of 3.5 
Mpc to be consistent with the dynamical studies mentioned above. 



2 OBSERVATIONS AND DATA REDUCTION 

We used the Ca mbridge Infra-R ed Panoramic Survey Spectro- 
graph, CIRPASS dParrv et ai]|2004l) . mounted at the Gemini South 
telescope, during the March 2003 queue scheduled observing cam- 
paign. Some general observational details are given in TableQ] 



Date 


Band 


Exp. time 


Seeing 


13. 03. 2003. 


J 


3x900 s 


Of/6 


19. 03. 2003. 


.1 


6x900 s 


ff.'S 



Notes : We adopt the date at the start of the Chilean night, rather than the 
Gemini default which would move all dates one day later. 

The IFU consists of a 490 element macrolenj]] array arranged 
in a broadly rectangular structure, with hexagonal lens packing. 
Through the use of re-imaging optics, in the focal plane, two lens 
scales where available at Gemini, f/.'36 and O'.^S, giving fields of 
view of 13 x4'.'7 and 9'.'3 x 3'.'5 respectively. We have used the 0725 
scale to study the central region of Cen A with the spectral res- 
olution of ~ 2.2 A/pixel (FWHM) and resolution element of 1.8 
pixels, in the spectral range between 1.1 — 1.35 /im. 

Data is recorded as series of on- and off-target exposures (each 
900 s). Off- target sky exposures were taken with a sufficient offset 
(600") from the Cen A nebulosity to avoid contamination from the 
galaxy. Data reduction was performed following the recipe outlined 
in the CIRPASS data reduction cookboofl The CIRPASS IRAF data 
reduction tasks where used along side a suite of custom software 
written in the IDL programing environment. 

Spectra are extracted form the beam-switched science frames, 
the sky frame and dome flats taken immediately prior to or after 
the science frames. Standard star frames are also extracted, with a 
beam switch pair placing the star at either end of the IFU to allow 
sky subtraction without the need for offset sky frames. An optimal 
extraction routine is used to allow for the crosstalk between fibres. 

Extracted spectra are then flat-fielded, using the normalised 
extracted flat field image. This step cannot be performed before the 
extraction as it would remove much of the information on the fibre 
profiles which we wish to extract. A wavelength solution and 2D 
transformation to a common wavelength scale is performed using 
the IRAF tasks identify, reidentify, fitcoo and transform. An average 
OH night sky spectrum is used for the wavelength solution, u sing 
only stronger lines from the line list of Maihar a et al.l d 19931) and 
free from blends. The average flat field response is then multiplied 
back into the data to suppress the spectrum of the flat field lamp. 

A relative flux calibration is derived using observations of the 
standard star Hip066842. For the star, a representative black-body 
spectrum is assumed and used to generate a spectral response func- 
tion and telluric absorption correction, after interpolating across 
strong stellar absorption features. The spectrum of Vega from the 
IRAF SYNPHOT packages is used as a reference for the magni- 
tude scale. Both the observed and the reference Vega spectrum are 
integrated over a clean region (away from telluric and strong stel- 
lar features) of the observed spectral range to find the calibration 
zero-point. 



2.1 CIRPASS data reduction 

The CIRPASS spectrograph, housed in an insulated cold room and 
cooled to -42° C, is a fibre fed instrument with observing capa- 
bilities in the J and H bands. Operating at a moderate resolution 
(R~3,000), the night sky OH emission spectrum is sufficiently dis- 
persed to allow one to capitalise on its intrinsic darkness between 

emission lines. 

CIRPASS has both integral field (e.g. ISharp et al1l2004l) an d 
multi-object observing modes (e.g. iDohertv et al. 20041 20061) . 



2.2 Mosaicing and alignment 

During typical CIRPASS observations, dither offsets on target are 
performed using integer lens offsets to allow straight-forward com- 



1 CIRPASS employs a macro lens array rather tha n a micro lens ar- 
ray, to avoid Focal Ra tio Degradation (FRD) losses l Ennico et af] 1 19981 ; 
iKenworthv et all 19981) 

2 The CIRPASS data reduction soft- 
ware and documentation is available from 
|http : / /www, ast . cam, ac . uk/~opt ics/cirpass /index .html| 
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bination of the observational data sets. However, data in each pass- 
band was taken using a number of different offset strategies. Align- 
ment and combination is complicated further by the hexagonal 
lenslet array. In order to align and stack the various data sets we 
therefore resample the data onto a regular Cartesian grid, deter- 
mine optimal image alignments via cross-correlation, and combine 
the data using integer pixel shifts. The process is as follows. 

For data taken on different nights the wavelength solution was 
significantly different. We therefore first correct all data in a spe- 
cific band to a common wavelength solution. We achieve this using 
an IDL implementation for the STARLINK FIGARO data reduction 
task SCRUNCH which allows accurate variance array propagation 
on rebinning. 

A Cartesian grid is generated which completely encloses the 
individual CIRPASS IFU frames. Since the entire Cen A data set 
was recorded using a position angle (PA) of zero, no account need 
be taken of rotation of the IFU coordinate frame. The pitch of the 
regular grid is chosen to oversample the IFU lens scale. Oversam- 
pling by a factor of three was found to give an excellent compro- 
mise between preserving the data integrity while keeping the output 
data set (which has of the order of 1000 wavelength channels) to a 
tractable size. A spectrum is assigned to each element of the regular 
array based on proximity to the elements of the IFU array. 

A set of integrated images is generated from these data 
cubes by summing along the dispersion axis, avoiding wavelength 
ranges subject to questionable data due to spectral curvature. These 
summed images are t hen cross-correla ted, using elements of the 
IDL astrolib libraries jLandsmarj[l993l) . Initial shift estimates are 
generated by eye from the summed images. The generated shifts 
are then used to apply integer pixel shifts to the resampled data 
cubes. 

Variation in the recorded intensity of spectra within the in- 
dividual observations was at the 10% level during each night and 
higher between nights (due primarily to the difficulty associated 
with observing accurate standard stars, over a range of airmasses, 
before, during and after each data set is recorded). We therefore 
derive a multiplicative scaling factor to apply to each data set be- 
fore the data is combined. The scaling factors are derived through 
measurement of the mean quotient between each data set and a ref- 
erence data set, for all spectra which are found in the complete 
overlap region between all IFU pointings (~33% of the data). 

After scaling each data cube, the data are mean combined us- 
ing a 3-sigma clipping algorithm to reject cosmic rays, which are 
still present to some extent in the data at this stage (as there are 
too few degenerate frames for a simple sigma clipping of the unex- 
tracted data) both as positive and negative features in the data due 
to the beam-switched sky subtraction. The resulting pixels scale is 
a factor of 3 smaller than the actual spatial resolution. We return 
to the original CIRPASS pixel size of 0'.'252 by averaging spectra 
within neighbouring 3x3 pixels. The error arrays were propagated 
by averaging in quadrature. 



signal Si and a measurement uncertainty (noise) Ni, then the S/N 
of the spectrum of n pixels is given by 



(S/N) = 



1 v-m o 



(i) 



In this approach the main contribution to the signal comes from the 
continuum (stellar) light and it is used for measurement of stellar 
kinematics. 

On the other hand, emission-lines do not necessarily follow 
the level of the stellar continuum and when estimating their S/N 
it might be more useful to estimate the 'signal' of the emission- 
line itself compared to the local noise of the spectral region sur- 
rounding the line. In this way, one is not including the information 
on the continuum and that is reflected in the res ulting bins, which 
then direct l y trace the emission-lines. Following iFathil d2005h and 
ISarzi et al] ( feOOtj) . we define amplitude-to-noise ratio (A/N) for 
each emission-line (continuum subtracted spectra), where the am- 
plitude (A) of an emission line is estimated as a median of the three 
largest pixel values in the emission-line, and the noise (N) is de- 
fined as the standard deviation of the spectra in the spectral region 
not contaminated by the presence of emission-lines. 

Since emission-lines are not necessarily present everywhere 
over the field of view, we used both binning approaches. The S/N 
based binning is used to present the data, analyse the global proper- 
ties and compare them with other observations (Section[3}. We used 
a target S/N of 20. The A/N based binning is used to measure the 
kinematics of individual lines (Section [4). We exp erimented with 
different ratios and found, similarly to dSarzi et alj 2006), that A/N 
= 2.5 is the most appropriate to detect lines and reduce the loss of 
spatial resolution. Both S/N and A/N were chosen such as to min- 
imise the degradation of spatial resolution in the central ~ 1". 



2.4 Additional data 

In addition to CIRPASS data, we also used other data sets to sup- 
plement our analysis. These include ground-based optical and ra- 
dio observations and space-based optical and X-ray imaging. We 
used reduced images of Cen A obtained by the Very Large Array 
(VLA) in A-array configuration and with the Chandra X-ray Ob- 
servatory, kindly provided by M. H ardcastle. Both images were 
observed, reduced and presented in iHardcastle et all d2003l) . We 
also used archival near-infrared observations of the central region 
of Cen A with the Hubble Space Telescope (HST). These include 
broad and narrow band imaging using NICMOS camera 2 (0'.'075 
pixel -1 ) with F187N (Paschen a) and F2 22W (K-band) filters . 
Th ese images were prev iously presented in ISchreier et ail (l998) 
and Marconi et al. I d2OO0t) . Finally, in order to constrain the mass at 
a larger radius we used a 2MASS K-band image of the galaxy. 



3 RESULTS 



2.3 Binning 

A minimum signal-to-noise ratio (S/N) is required to reliably mea- 
sure the properties of spectral lines. To achieve this, we spa- 
tially binned our merged data using an adaptive algorithm of 
ICappellari & CopirJ (2003) creating compact and non-overlapping 
bins with uniform S/N in faint regions. A desirable value of S/N de- 
pends on the purpose of the analysis. Similarly, a way to estimate 
the S/N can be different. Usually, if each pixel in a spectrum has a 



Our observations cover the nucleus of Cen A. The orientation and 
the size of the merged CIRPASS field can be seen in Fig.Q] In this 
section we present the CIRPASS data cube and compare it with 
other observations. 



3.1 Nuclear spectrum 

All spectra within a circular aperture of 0'.'5 radius were combined 
to form the nuclear spectrum presented on Fig. [2] There one can 
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Table 2. Emission-line ratios within circular apertures of different sizes cen- 
tred on the nucleus 
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Figure 1. Broad band observations of central region of Cen A. Left: ground- 
based 2MASS K-band image; the over-plotted square shows the size and 
orientation of the right image. Right: HST/NICMOS F222M image with 
over-plotted footprint of the CIRPASS field-of-view. On both images North 
is up and East to the left. HST image has been rotated accordingly. 
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Figure 2. Total spectrum from a circular aperture of 0'.'5 radius, centred 
on the nucleus. Three emission-lines used in this study are labelled at the 
wavelength redshifted for 543 km s — 1 . The vertical dashed lines show the 
position of two candidate lines contributing to the blue wings of detected 
emission-lines. There is no evidence for existence of HeIA11970 line (dot- 
dashed line). The vertical lines at the bottom of the spectrum show positions 
of the sky lines in the observed wavelength region. 



see the three major emission-lines which are of interest for this 
study. They are, in order of their observed relati ve strengths: [Fell l 
A12567, Paschen (3 A12818 and [PII] A11882 dOliva et : alj[l99(ih . 
The observed lines are redshifted, with respect to their laboratory 
values, by 22 A, corresponding to a vel ocity of 543 kms" 1 . The 
first two lines were already observed bv lMarconi et al. I d200lh hav- 
ing the ratio of ~ 1.7, but the observation of the [PII] line is, to our 
knowledge, the second extragalactic measuremen t of thi s emission- 
line, the first being in NGC 1068 by lOliva et ail d200ll) . The blue 
wings of the [Fell] and Pa/3 emission-lines could be attributed to 
Hel A12530 and Hel A12790 emission. There is, however, no evi- 
dence for presence of HeIA11970 emission, which suggests Hel is 
generally not strong. The origin of the red wing of the [PII] is not 
known. The detected emission-lines are not contaminated by sky 
lines, except a minor contribution from the sky lines to the bluest 
part of the Pa/3 wing. 

The relative strengths of the lines can be used to determine 
the nature of the ionisation mechanism. We estimated the strength 
of the emission by fitting double Gaussians (see Section |4~T1 for 



aperture 


[FeII]/Pa/3 


[FeII]/[PII] 


(1) 


(2) 


(3) 


0.25 


2.3 


3.0 


0.5 


2.8 


3.5 


0.75 


2.1 


3.7 


1.0 


2.2 


1.2 


1.25 


2.6 


4.7 



Notes - Col.(l): Radius of aperture [arcsec]; Col. (2): ratio of [FeII]/Pa/3 
flux Col.(3): ratio of [FeII]/[PII] flux. 

more details) within apertures of increasing diameter, centred on 
the nucleus. Table[2]shows the measurement of the ratio of the lines. 
Within a 0//25 aperture [FeII]/Pa/3 ~ 2.3, while [FeII]/[PII] ~ 3.0. 
As the aperture increases the [FeII]/Pa/3 ratio fluctuates between 
2. 1 and 2.8, being typical of active galaxies or supernova remnants, 
as c an be seen from observations and line diagnos t ic diagrams 
(e.g. [Simpson etalJI 19961: lAlonso-Herrero et al ] |l997l : lMourietal.l 
l200fAlRodriguez-Ardila et alj2004l) . The active nucleus in Cen A is 
an obvious source of [Fell] emission-lines, and although it is possi- 
ble that a fraction of the [FeII]/Pa/3 ratio is caused by star-formation 
(e.g. supernova remnants in the star forming regions), the fact that 
[Fell] emission is associated with shock, also points towards a mass 
outflow from the nu cleus (e.g. the jet) a nd resulting shocks as the 
source of ionisation dKawara et al.ll 19881) . 

The [FeII]/[PII] ratio, however, shows a clearly increasing 
trend, reaching 4.7 at the largest aperture. This can be used to asses 
the im portance of shock s in the total ionisation mechanism of [Fell] 
lines. lOliva et al.l d200ll) found that the [FeII]/[PII] ratio is a good 
indicator for distinguishing the origin of the emission-lines, espe- 
cially if they are ionised through photoionisation by soft X-rays or 
shocks propagating through the interstellar medium. The authors 
state that in regions with [FeII]/[PII] < 2 gas is photoionised, while 
if [FeII]/[PII] > 20 the material is excited by collisions in fast 
shock waves. The observed values of the ratio are between these 
two limiting values, although closer to the upper limit for pho- 
toionisation, suggesting the AGN X-rays emission as the source 
of photoionisation. Still, the rising values argue that shocks play an 
increasingly important role in ionising the interstellar medium fur- 
ther away from the nucleus. This could be reflected in the motion 
of gas clouds, where a more settled motion should be observed in 
the central regions. 



3.2 Extent of the emission 

Figure [3] presents the spatial extent of the I-band CIRPASS obser- 
vations of Cen A. We binned the final merged data cube to an S/N of 
20. The binned spectra were continuum subtracted after fitting 4th 
order polynomials to emission-free parts of spectra. The emission- 
lines presented in Fig. [3] were measured by summing the contin- 
uum subtracted spectra along the wavelength dimension in limited 
spectral regions centred on the lines themselves (see the caption of 
Fig.[3]caption for details). 

Maps show two spatially separated regions of high integrated 
flux, one in the nucleus (for the determination of the central pixel 
see Section [3~3l l and the other in the North-East region (NE), about 
2"away from the nucleus. The detected emission-lines are mostly 
confined to the central 1", but Pa/3 and [Fell] emission also extend 
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Figure 3. Maps obtained by summing the data cube along the spectral di- 
mension after binning to S/N = 20. The maps are oriented such that North 
is up and East to the left. From left to right: summed spectra with over- 
plotted VLA radio observations, spectra summed between 12554 - 12622 
A ([Fell] line), spectra summed between 12807 - 12874 A (Pa /3 line) with 
over-plotted Paa emission from archival HST/NICOMS observations, spec- 
tra summed between 1 1870 - 1 1937 A ([PII] line). Note the change in the 
scale for the [PII] map, the agreement between Paa and Pa/3 structures and 
the coincidence of the Nl knot in the jet with the North-East flux peak in 
our observations. Flux F x is in units of 10 — 17 ergs/s/m 2 /A. 



towards the NE. Pa/3 is also present in two blobs about 3"away 
from the nucleus. The weakest emission comes from [PII]. This 
emission-line falls off rapidly and outside 1" the summed spectra 
are dominated by noise, but one can obtain information on some- 
what larger radius (e.g. within 1'.'25) by fitting the spectral line as 
done in Sectionl4~T1 



3.3 Comparison with observations on different wavelengths 

In order to determine the position of the nucleus in the CIR- 
PASS data, we com pared our observations with HST imaging. 
ISchreier et al.l dl998h identified the nucleus with the peak in Paa 
emission, which also corresponds with the peak visible on K band 
(F222W) NICMOS images and the position of the radio nucleus. 
We find the position of Cen A's nucleus in our observations by as- 
suming that the Paa and Pa/3 emission have the same spatial origin 
and by cross-correlating the NICMOS F187N narrow band image 
with an image constructed from our unbinned data cube summed 
over the spectral region of the Pa/3 emission. 

The spatial distribution of the two Paschen lines matches well, 
as it can be seen in Fig. [3] where the NICMOS Paa im age is over- 
plotted on the binned Pa/3 image. ISchreier et al. determined 
the position angle of the elongation to be ~ 33°. Also, the two off 
centred blobs are well matched in this comparison. The extended 
distribution of Paa emission is also traced by Pa/3 emission. 

We also compared the CIRPASS obs ervation of Cen A wit h 
existing radio and X-ray observations dHardcastle et al. I l2003h . 
Most of the action in the radio and X-ray happens at larger radii. 
This is especially the case for X-rays where the central region is 
featureless, and the very centre is saturated. The existence of X- 
rays, however, suggests a possible excitation mechanism for the 
observed spectral lines. 

Radio observatio ns trace a jet with a position angle of 
51° jClarkeetalJll992l) down to the radio nucleus. Radio contours 
are over-plotted on a CIRPASS image obtained by summing all the 
spectra along the wavelength dimensions (Fig. [3}. The direction of 
the radio jet coincides in projection with the NE region of high 
flux. The jet crosses the region of high flux loosing its energy: at 
the start of the region the energy in the jet is 10 mjy/beam and at 
the end the emitted energy is about 1 mJy/beam (the lowest plotted 
contour). This drop in jet's intensity was used to define knot Nl by 
iBurns et all dl983l) . T he direction of the jet is about 18° away from 
the orientation of the Schreier et a L1( ll998h Paa disk, but between 
the two Pa/3 (and Paa) blobs, touching the top of the East blob at 
the position of the kink in the jet. The kink in the jet could be ex- 
plained as a deflection of the jet from a cloud complex which is 
now glowing in Pa/3 light. 



3.4 A hole in the extinction map or a jet-induced star 
formation region? 

The most striking feature of the reconstructed images is the sec- 
ondary peak in the integrated light which, in projection, is co- 
spatial with the radio jet, bounded by the two Pa/3 blobs and arises 
at the end of knot Nl. Cen A was the first object for which it was 
suggested that star formation in certain r egions was jet-induced 
dGraham & Pricdl 198ll: iBrodie et al.lll983l). soon followed by the 
Minkowski object dvan Breugel et alj|l985h . and numerous other 
objects with different environments and jet powers at low and high 
redshifts. The connection between the suggested nuclear star for- 
mation and the radio jet is, however, not as straight forward as at 
large scales, where blobs of optical emission and enhanced contin- 
uum correspond to specific jet structures. In this case, the end of a 
jet knot is followed by an increase in the total near-infrared flux, 
mostly coming from the continuum (in the same region there are 
no strong emission-lines), and bounded with regions of detectable 
emission: North and East Pa/3 blobs and the nucleus. Interestingly, 
the East Pa/3 blob coincides with the kink in the radio jet. 

There are a few possible explanations for these observations. 
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Figure 4. An example of two spectra, one from the North-East region of 
high continuum level (top spectrum) and the other from the opposite South- 
west region with typical a continuum level (bottom spectrum). The spatial 
coordinates are (—1", — 1'.'26) and (l'.'l, — 1'.'7), for top and bottom spec- 
trum, respectively. 



The two most likely are that the NE flux peak is related to a jet- 
induced star formation or that it coincides with the location of a 
hole in the extinction. In the first scenario, star formation could 
happen under the influence of the jet as it plows through the inter- 
stellar medium pushing gas clouds away from its path, and result- 
ing in a region with higher stellar continuum. In the same time, Pa/? 
blobs could be ionised from the central engine or the jet, provided 
that they are on the unobscured path of non-thermal photons. In 
the second scenario, the NE flux peak is visible because it is un- 
der an extinction hole. In addition to that, the North and East Pa/3 
blobs may also originate in star forming regions where young stars 
provide enough photons to ionise gas clouds. 

The extinction hole is supporte d by HST/WFPC2 F 814W im- 
ages (third panels on Figs. 5 and 8 of lMarconi et alj2OO0l) on which 
the nucleus (a faint point source) and the less obscured NE region, 
which looks like a low extinction hole, are separated by a dust lane. 
The dereddend F814W image does not have similar features. The 
F222M image (K-band, presented on the same figures), which is 
less influenced by obscuring material shows only a marginal in- 
crease of flux in the NE region. The spectra from NE (high contin- 
uum) and SW (low continuum) regions are, however, very similar 
(Fig. |4), NE spectra being systematically higher, without any sug- 
gestion of change in the spectral slope, as it would be expected from 
regions of different extinction. 

Assuming that jet induced star formation is the source of the 
high continuum and the Pa/3 emission in the blobs, it is possible 
that the North Pa/3 blob was also once induced by the jet (the jet is 
currently on the East Pa/3 blob), which precessed down with time 
(in projection) from PA =~ 25° to present PA =~ 50°. 

It is hard from the present observations to disentangle the true 
origin of the properties of the North-East region. The precession of 
the jet is perhaps a far-fetched assumption, but our data, which are 
not complete enough for a decisive conclusion, support marginally 
more the jet-induced star formation in the continuum region, where 
an extinction hole can also contribute to the detection. 



4 KINEMATICAL PROPERTIES OF THE NUCLEAR 
DISK 

In the previous section we binned to a S/N ratio of 20 to present 
the data and discus the main trends. In this section we concentrate 
on the three observed emission-lines and use the A/N ratio based 
binning. We used a value of A/N = 2.5 for all three lines. 

4.1 Narrow and broad components 

The shape of the emission-lines on Fig|2]suggest that single Gaus- 
sians are not always adequate line models. Hence, we fitted the 
spectral features with double Gaussians when necessary. Each of 
the three emission-lines in a given spectrum were fitted separately 
and independently. Firstly, we fit the spectra with a 4th order poly- 
nomial avoiding the regions with the emission-lines. We subtract 
this continuum contribution and each emission-line we fit with a 
single Gaussian. These two fits are used as preparation steps for 
the final fit of the emission lines. In the final fit, we fit a double 
Gaussian and a 4th order polynomial simultaneously to the original 
spectra (not continuum subtracted). The parameters of the contin- 
uum fit and the single Gaussian fit are used as initial values in the 
minimisation procedure of the final fit. The results of the contin- 
uum fit were taken as the initial values of the 4th order polynomial 
parameters. The initial values of the double Gaussians parameters 
were chosen from the single Gaussian fit in the following way: (i) 
sigma (width) of the first Gaussian was set to be a half of the sigma 
of the single Gaussian fit; the sigma of the second Gaussian was 
kept the same as in the fit with the single Gaussian, and (ii) the 
mean of the first Gaussian was set to be equal to the mean obtained 
by the single Gaussian fit. The mean of the second Gaussian was 
set to be about 730 km/s (30 A) less than the systemic velocity of 
the galaxy. 

To robustly constrain the final fit we introduced the following 
limits for the parameters of the first Gaussian: (i) the mean was al- 
lowed to vary 250 km/s around the initial value, and (ii) the sigma 
was bounded between 4 and 15 A(or 95 km/s and 350 km/s). The 
lower boundary is the instrumental resolution. The second Gaus- 
sian had the following constraints: (i) the mean was allowed to vary 
500 km/s around the initial value, and (ii) the sigma was bounded 
between 4 and 45 A(or 95 km/s and 1000 km/s). The upper bound- 
ary was seldomly reached. The final fit with a do uble Gaussian was 
performed using the MINPACK implementation ( Mor e et aljl98d) 
of the Levenberg-Marquardt least-square minimisation methoqj. 
The double Gaussian system increased the goodness of fit, but it 
was not required at all positions for all emission-lines. The sec- 
ond Gaussian was needed to fit the Pa/3 emission-line almost ev- 
erywhere in the field while for [Fell] emission, it was necessary 
only in the centre, and for [PII] almost not at all. 

Errors on the parameters of the Gaussians (kinematics) were 
obtained through Monte Carlo simulations. For each double Gaus- 
sian fit to the data, we derived uncertainties from 100 random real- 
isations of the best fit, where the value at each pixel is taken from a 
Gaussian distribution with the mean of the best fit and standard de- 
viation measured on the residuals between the fit and the data. All 
realisations together provide a distribution of values from which we 
estimate la (68.3%) uncertainties. 

It is not necessary to expect that binning will always reach the 



3 We used an IDL version of the code written by Craig B. Markwardt and 
available from: http://astrog.physics.wisc.edu/~craigm/idl 



IFS of Centaurus A 7 



Pa f [Fell] 




-3-2-10 1 2-3-2-10 1 2 



Figure 5. Maps of total intensity and kinematics of detected (A/N = 3) 
emission-lines. Left: Pa/3; right: [Fell]. From top to bottom: the total flux 
of the emission, the mean velocity (centre of the Gaussian), the velocity 
dispersion (<r of the Gaussian, without the contribution from the instrument 
dispersion). Over-plotted contours are the contours of the total flux. Dif- 
ferent symbols relate to spectra in Fig. [6] Flux F x is in units of 10 — 17 
ergs/s/m 2 /A. Regions without detected emission are shown in dark grey. 

set levels, since the emission need not to be present, and binning the 
noise will not create signal. Calculated amplitudes of both narrow 
and broad Gaussians and the standard deviation of the residuals of 
the fit can be used a posteriori to determine the true A/N ratio 
of our data. Results show that the reached A/N drops below the 
targeted level outside of the central region (size depends on the 
emission line). The quickest spatial drop in A/N is seen for [PII], 
and the slowest for [Fell] emission-lines. Also, the highest A/N is 
meas ured for [Fell] an d the lowest in [PII] lines. 

ISarzi et all |2006) extensively discuss at which level of A/N 
ratio an emission-line is robustly detected. This depends on the lo- 
cal spectral features (existence of complex absorption features will 
require a higher A/N ratio), but, generally, an emission-line with 
A/N = 2.5 — 3 can be considered as detected. We binned to 
A/N = 2.5 to keep the spatial resolution as high as possible, but in 
a posteriori check we chose a more conservative limit of A/N = 3, 
and in the rest of the analysis we consider only emission-lines 
above this level. 

Emission-lines that met the above criterion are presented in 
Figs.[5]and[6] The two rows of Fig.[5]show the two-dimensional ex- 
tent, morphology and kinematics of Pa/3 and [Fell] emission-lines, 
respectively from top to bottom. The [PII] emission-line reached 



A/N = 3 level only in the central (/.'5 and the quality of the ex- 
tracted kinematics was rather poor, so we do not consider it in the 
rest of the study. Fig.|6]shows examples of the spectral features ob- 
tained from the same spatial positions indicated by the symbols on 
both figures, but for different emission-lines (from top to bottom: 
Pa/3 and [Fell], respectively). 

Although at almost all positions a combination of a narrow and 
a broad Gaussian is required to fit the Pa/3 emission-line, only in a 
few cases the broad component can be really recognised as a sepa- 
rate line. More often, the broad component serves as an additional 
continuum correction (Fig. |6). When present, the broad Gaussians 
are blue shifted with respect to the narrow Gaussians for (alway fl) 
795 kms" 1 . The spatial distribution of the broad component (with 
A/N > 3) is limited and concentrated in the very centre and some- 
what elongated in the NE direction, near the position of the jet and 
between the blobs of the narrow component. Line widths are gener- 
ally less than 800 kms -1 (sigma of the Gaussian), and as such are 
not comparable to the emission from classical broad-line regions 
in AGNs. We call this component broad only to distinguish it from 
the main line emission and we do not consider it further. The main 
(narrow) component is detected on a larger area in the centre and at 
the position of the Paa blobs, (Fig. [3j. Spatially, the distribution of 
the narrow Gaussians seems to be elongated along the NE direction 
with PA ~ 30° . Thi s value is consistent with the position angle of 
ISchreier et all i 19981) Paa disk. 

Different results are recovered from the [Fell] line. The broad 
Gaussian component to the [Fell] emission-line is almost com- 
pletely absent, as expected for a forbidden line, except in the very 
centre where A/N is still below 3 and broad Gaussians improve the 
fit to the continuum more than fit a real (robust) emission feature 
(Fig.[6j. The [Fell] line is detectable on a larger area than Pa/3 and 
it is not elongated, but more circularly distributed. The emission is 
concentrated in the centre, but also spread towards the NE corner, 
approximately following the position of the jet, although it is just 
above the detection limit and patchy distributed, mostly over the 
East Pa/3 blob and approaching the position of the North blob. 

4.2 Kinematical evidence for a disk 

Both Pa/3 and [Fell] emission-lines show rotation and they are 
likely part of the nuclear disk structure discovered bv lSchreier et al.l 
dl998h . Here, we wish to describe their kinematic properties as they 
are given on Fig. [5] 

The there regions with measured Pa/3 emission (the centre, the 
East and the North blob), although spatially separated, seem to be 
part of the same kinematic structure. The central blob (covering 
roughly the central 1 arcsec 2 ) shows a quite regular rotation, with 
zero-velocity curve (ZVC) lying horizontally in the East - West 
(EW) direction. The East blob is predominantly blue-shifted, while 
the North blob is red-shifted, with respect to the systemic velocity 
(the median velocity across the map is ~ 560 kms" 1 ). All three 
blobs have similar velocity amplitudes and there is a clear division 
on the map: North is receding and the South(-East) is approach- 
ing the observer. If the emission were detected at all positions, it is 
likely that we could see a continuous kinematic map, with a twist- 
ing ZVC: from horizontal (EW direction) to almost 45° (NE-SW 
direction). 

Velocity dispersions, measured as the width (a) of the narrow 

4 We also fitted the line by fixing the relative distance of the Gaussian 
peaks, which gave equal results with marginally worse \ 2 - 
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Figure 6. Example of spectral features recovered from different positions in the nucleus of Cen A for different emission-lines. Top: Pa/3; bottom: [Fell]. 
Symbols in top left corner are the same as symbols in Fig. \5\ Every symbol shows the spatial location of the bin from which the spectrum was taken. Only 
regions around the given emission-line used in the fit are shown. Yellow, blue and light blue dashed lines are narrow, broad Gaussian components and the fit to 
the continuum, respectively. The red line shows the combined fit to the data. Note the strong contribution of the broad Gaussian to the central Pa/3 emission. 
Error bars (la) are plotted for every 10th data point for clarity. 



component Gaussian, confirm the relations between the blobs. The 
East and North blob have low velocity dispersion as well as the 
outer regions of the central blob. The dispersion increases towards 
the centre (staying below 200 kms -1 ), possibly being elongated 
perpendicular to the orientation of the central blob (PA ~ —60°). 
The low velocity dispersion over most of the region also suggests a 
regularly rotating structure, such as a thin disk of gas clouds. 

[Fell] emission covers a larger area and it fills some of the 
gaps seen on the Pa/3 maps. The distribution of velocities is remark- 
ably similar between the two lines. [Fell] emission also shows rota- 
tion (the median velocity across the map is ~ 590 kms -1 ), which 
changes direction: the ZVC is oriented E-W in the centre, but be- 
yond 1" it twists towards NE-SW. Actually, the ZVCs on Pa/3 and 
[Fell] maps have identical spatial positions and orientations, where 
the North on [Fell] map is receding and the South is approaching. 
Beyond 2" where the emission becomes patchy, it is still possible 
to trace the rotation of the whole structure. 

The velocity dispersion is low in the outer regions, but rises 
towards the centre. It is higher than the Pa/3 velocity dispersion 
everywhere, except at the edges of the detection region. It also rises 
faster (median level is ~ 127 kms - ), until levelling at a plateau 
(~ 200 kms" 1 ) which is elongated in the same directions as the 
enhanced Pa/3 velocity dispersion structure. The orientation of the 
plateau (PA = —60°) is unexpected also because it is not aligned 
with the the ZVC. An increase along the ZVC could come from 
unresolved rotation in the steepest part of the velocity map (where 
velocities change from approaching to receding), but even in this 
extreme case the velocity dispersion would not be as elongated. We 
will follow this issue further in Sections [4.3l and l531 

From Fig.[5]it is obvious that different blobs of Pa/3 and [Fell] 
emission have related kinematics, even when they are spatially 
separated. A similar finding is also visible on SAURON maps of 



ionised gas in early-type galaxies i Sarzi et alj |2006h and bulges 
of Sa galaxies dFalcon-Barroso et aL 20061) , where gas although 
patchy or along filaments shows continuous kinematics, suggesting 
that all emission elements belong to the same dynamical structure, 
but that the level of ionisation, the amount of gas or extinction is 
not the same and varies locally. 

The rotation of both lines is to a first order regular and sug- 
gests that the emission is distributed in a disk-like structure. It is 
also clear that the kinematics of Pa/3 and [Fell] emission are gener- 
ally similar (velocity orientation and amplitude), but there are some 
notable differences such as a higher velocity dispersion in [Fell] 
lines. Judging from velocity dispersion maps, it seems that [Fell] 
emission is more perturbed than Pa/3 and its motion is not com- 
pletely ordered nor purely disk-like. Deviations could be caused 
by physical processes such as inflows, outflows, non-axisymmetric 
perturbations in the potential, warps or simply material not belong- 
ing to the disk. We investigate this in more details below. 

4.3 Orientation of the velocity structures 

Assuming that the velocities represent material moving in a thin 
disk, it is, in principle, possible to determine both PA and incli- 
nation (i) at the same time and as a function of radius. This is, 
in practice, often quite difficult, especially if there are deviations 
from circular rotation or data are noisy. In our case, we decided to 
determine only global properties and to do that through a two stage 
process. 

We first follow the method for determination of th e global 
kinematic PA from Appendix C of lKrainovic et al. Briefly, 
one constructs a set of bi-anti-symmetric velocity maps (because 
of the assumed axisymmetry) with different PAs and compare 
them with the observed velocity map. The likeness of a model 
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Figure 7. Maps of the observed mean velocity (far left), circular velocity (middle left), residuals between the mean and circular velocity (middle right), and the 
best fitting thin disk models (far right) using the derived orientation for Pa/3 (top) and [Fell] (bottom) emission-lines disks. Circles over-plotted on observed 
velocity maps have radii of 0'.'9 (blue) and l'.'S (black). Bins marked with black dots were not included when constructing circular velocity maps. The plus 
symbol denotes the origin of the map, assumed to be the position of the Cen A nucleus. Standard deviations of residual maps are 18 km s — 1 and 25 kms 
(but 16 kms -1 within the blue circle), for Pa/3 and [Fell] maps, respectively. The method to construct circular velocity maps is described in Section l43l and 
the thin disk model is described in Section[5] 



map (V m ) to the observed map (V) is measured by the quantity 
X 2 = Tl^ \{Vm,j — Vj) / AVj] 2 , where AVj are measurement er- 
rors on velocities for each bin j on a map m. The model map with 
the PA closest to the observed value will contribute with the small- 
est x 2 - 

We applied this method to [Fell] and Pa/3 maps which were 
slightly modified by selecting only spatially continuous emission 
on the maps. This meant decreasing the size of maps to the central 
regions (about 1 arcsec 2 for Pa/3 and about 2 arcsec 2 for [Fell]), but 
increased the robustness of the method. For these trimmed maps 
(see the first column in Fig. [7j we obtained the following values: 
PA = -3° ± 10° for Pa/3 and PA = -23° ±3° for [Fell] velocity 
maps. All errors quoted here are at 3u confidence level. The [Fell] 
map is about 2 times larger than the Pa/3 map and difference in PA 
comes from the change in the ZVC as one includes a larger area. 
On Fig. |7] the blue and black circles on the maps show the regions 
used in the fit to determine the PA for Pa/3 and [Fell], respectively. 
Repeating the exercise for [Fell] within the blue circle, one gets 
PA = —14° ± 5°, which is in good agreement with the Pa/3 value. 

The next step is to estimate the inclination of the disk assum- 
ing that the motion in the disk is purely circular. In that case, a 
velocity profile extracted along an ellipse will be described by a 
cosine law, if the ellipse samples equal radii in the disk plane and 
has the same PA as the disk. This ellipse is then related to the in- 
clination thr ough its elliptici t y: 1 — e = cos i. Using the kinemetry 
software of Krai no vie et alj ( I2006T) and assuming a constant PA 
measured above, we estimate the i of the maps in the following 
way. Starting from a given i we stepped through the velocity map 
extracting velocity profiles along ellipses with 1 — e = cos i and in- 
creasing semi-major axis 7?. Each velocity profile we least-squares 



fitted with the formula V(R, 8) = ( V C (R) / sin i) cos 6, where V c 
is the circular velocity and 6 eccentric anomaly, measured from the 
projected major axis of the ellipse (given by the fixed PA). We 
linearly interpolated V(R, 9) to V(x, y) creating a model velocity 
map at the observed coordinates on the sky. This has been done for 
a range of different inclinations and each model velocity map (V m ) 
was compared to the observed map (V) similarly as above using a 
quantity: x 2 = ^2f[(V m ,j — Vj)/^Vj] 2 , where AV are measure- 
ment errors. The inclination of the model map with the minimum 
X 2 was recognised as the global estimate of the inclination of the 
observed disk. An almost identical procedure (they varied both P A 
and i simultaneously) was used previously bv lBuson et alj f2006). 

In the case of Pa/3 emission, setting the PA to —3° we ob- 
tained i = 25° ± 5° . Setting the PA to -23° for [Fell] disk we ob- 
tained i — 30° ± 5° . The two inclinations a re consistent w i th eac h 
other, and a re als o consistent with studi es [Marc oni et all ( 1200 lh . 
ISilge et alj feOOj^ iMarconi et al. (2006)), but somewhat smaller 
than models by JHaring-Neumaver et alj|2006l) . In the second col- 
umn of Fig. |7] we show the reconstructed circular velocity maps 
using the best fitting PA and i. 

One can most clearly see the deviations from a pure circu- 
lar motion on the residual maps obtained by subtracting the recon- 
structed circular velocity maps from the observed velocities (shown 
in column 3 on Fig. [Vj. Both residual maps have visible non- 
zero components. On average, residuals are higher on the [Fell] 
map (standard deviation is 18 kms -1 for Pa/3 and 25 kms -1 for 
[Fell]), but within the central region (blue circle) the residual are 
very similar and have the same spatial trend )standard deviation of 
~ 16 kms -1 ). This is consistent with the rising [FeII]/[PII] ratio 
towards the edge of the map, where shocks become an increasingly 
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Figure 8. Measured emission-line velocities and model predictions (see 
Section \5\ for details) along PA = 0° . The [Fell] velocity data and the 
best fitting model are shown with triangles and the dashed-dotted line, re- 
spectively. The Pa/3 velocity data and the best fitting model are shown with 
filled circles and the black line. 



important factor in the ionisation mechanism, and possibly influ- 
ence the kinematics. Shocks can result from the jet interaction with 
the surrounding matter, but can also occur between the tilted rings 
of the warped disk. On larger scales ([Fell] residuals map) there is 
a hint of a more complicated pattern which primarily arises from 
the difference between the observed (changing) and assumed (con- 
stant) orientation. This, possibly three-fold, pattern can be associ- 
ated with inflows, outflows or warps. 

The most distinct feature on the velocity maps is the twist of 
the ZVC bey ond 1". Inside this r adii the black hole is dominating 
the potential (M arconi et alj2 006) and is able to stabilise gas clouds 
in a disk. The central blob of Pa/3 emission is largely confined to 
this region only, but iron is ionised in almost continuous gas clouds 
up to 2", and the influence of instabilities is visible in the central 
emission-line disk. The two other Pa/3 blobs at ~ 2.5" — 3" from 
the nuclei, behave similarly like [Fell] disk on large radii, being far 
from the sphere of influence of the black hole where the instability 
can perturb the disk. 

Th e dusty disk in C e n A is known to b e warped at kpc 
scales jBlandetalj Il987l ; lOuillen et ail 1 19921; Nicholson etail 
1992; Ouill en et alj 1 19930. There are evidences for a bar spiral 
system ^Mirabel et al]|l999l : iBlock & SauvageH200Cl : iLeeuw et al.1 
2002), although the recent infrared observations support the warp 
disk hypothesis dOuillen et alj2006» . Our data shows that the warp 
(or a bar-like perturbation) can be traced all the way to 20 pc scales. 
It is, however, difficult to be conclusive on the matter from these 
data, and the meaning of residuals has to be taken with care, es- 
pecially when the measured uncertainties are taken into account: 



and lOuillen et al.1 ( |2006|) . we suggest that the origin of the high ve- 
locity dispersion is in the unresolved rotation coming from differ- 
ent folds of the warped disk. We draw this conclusion from Figs. 10 
and 11 of lNicholson et al]dl992l) and Figs. 5 and 7 of lOuillen et al.l 
(2006), which show that the disk changes its orientation such that 
in the inner part of the disk (r < 60") its SW points and in the outer 
part (r > 100") its NE points are closest to the observer. Although 
both studies construct the warped disk model to larger radii (beyond 
40") than what we present in this work, the orientation of the high 
velocity dispersion plateau (PA = —60° or PA — 120°) agrees 
well with the measured orientation of the warped disk. This means 
that our line-of-sight crosses the disk a few (at least two) times, and 
the observed kinematics will depend on where the ionisation of the 
gas clouds occurs. 

The [Fe II] emission-line arises from partially ionised regions, 
i.e. regions with a few per cent fraction of HI in the HII state. These 
regions are typically seen in the narrow line region environment of 
an AGN, where X-ray photons from the central engine can pen- 
etrate deep into the cloud. Thus, partially ionised regions do not 
have a direct view of the central engine, being shielded by a col- 
umn of ~ 10 20 cm -2 of Nh- The observed [Fell] emission could 
originate in different folds of the warped disk, also further away 
from the central engine. The resulting large difference in velocities 
along the line-of-sight, observed as unresolved rotation, has two ef- 
fects: the velocity gradient is shallower and the velocity dispersion 
higher than in the case of resolved rotation. 

In Fig. [8] we compare the central velocity gradient between 
the Pa/3 and [Fell] emission-lines, extracted from velocity maps 
along a slit of 3" length, C/.'25 width and PA = 0°. Within Of! 5, 
Pa/3 emission shows a steeper rotation curve than [Fell] emission, 
while around 1" both emission-lines reach the same level. In the 
substantially warped gas disk, the innermost regions would provide 
the required shielding for the [Fell] emitting material slightly fur- 
ther out. The lack of [Fell] emission from the innermost region, in 
contrast to Pa/3, would explain a slower rise of the [Fe II] line-of- 
sight velocity curve with increasing radius, plotted in Fig. [8] Also, 
unlike Pa/3, the [Fell] emitting material is confined to specific re- 
gions of the gas disk, increasing the likelihood that it provides an 
incomplete picture of the gas kinematics in the central region. The 
higher velocity dispersion indicate a more turbulent medium, one 
in which shocks may be playing an important role. This too, raises 
doubts about the [Fell] emitting gas serving as a good tracer of the 
gravitational potential in the innermost nuclear regions. 



16 km s~ 1 for Pa/3 and 19 km s~ 1 for [Fell] . We will, however, fol- 
lowing the recent results, continue in the framework of the warped 
disk. 

Comparing the orientation and the ex tent of the velocity dis- 
persion plateau with the tilted ring model o f lNicholson etaj]<1992l) 



5 EVIDENCE FOR THE CENTRAL BLACK HOLE 

In the previous section we showed evidence that the observed 
ionised gas is moving in a disk-like structure, which is consistent 
with being warped beyond 1". In this region, the motion is not com- 
pletely ordered, but the level of the residuals is comparable with 
measurement uncertainties. We also showed that [Fell] emission- 
line kinematics suffers from observed unresolved rotation, which 
can be explained with warped disk morphology. We continue the 
analysis by constructing simple dynamical models, based on the 
observed properties of the disk (orientation, circular motion) with 
a purpose of testing the hypothesis that there is a central black hole 
in Cen A that can explain the observed kinematics. 
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Table 3. The MGE parameters of the circular PSF. 
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0.041 
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0.126 
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0.080 


0.508 



5.1 Mass models 

As reported by Silge et al, (2005), the enclosed stellar mass within 
the central 1" is 3 x 10 7 M©, for a deredened profile, and 2 x 
1O 7 M0 for an observed (not deredened) light profile. The mass 
of the black hole exp ected from the M. — a relation (as given by 
iTremaine et alj200 2)is3x 10 7 M Q .The mass of the central black 
hole in Cen A determined by previous studies is, however, around 
1 x 10 8 Mp). As suming a velocity dispersion of ~ 150 kms -1 
dSilge et alfeOOl) . the sphere of influence of this object is ~ 1'. 3. 
This suggest that within 1" the dominant contribution to the poten- 
tial comes mostly from the central black hole. The total potential 
can then be written in the following way: 

$ = + (2) 

where the black hole contribution is simply given by $. = 
GM./r, while <E>* can be derived from the stellar luminosity den- 
sity assuming a mass-to-light ratio M/L. The stellar luminosity 
density of a system is obtained by deprojecting the observed sur- 
face brightness distribution, assuming a symmetry type (e.g. ax- 
isymmetry) for the system and its inclination. This can be con- 
venient ly done using the the Multi-Gaussian ex pansion method 
(MGE; lMonnet et al.Il992l ; lEmsellem et al.ll 1994b. 

Using the MGE software of ICappellaril d2002h , we simulta- 
neously fitted a ground-based K -band image from 2MASS Large 
Galaxy Atlas dJarrett et alj|2003l) ( 3"< r < 200") and the space- 
based NICMOS F222M image dSchreier et al J 1998b (0"< r < 9"). 
The 2MASS image was scaled to the NICMOS image, while the 
sky was measured on the 2MASS image. The MGE software ana- 
lytically deconvolves the instrumental PSF from the observed light, 
which we determined using the Tiny Tim software iKrist & Hookl 
d200ll) and also parameterised by three two-dimensional Gaus- 
sians, normalised such that t heir amplitudes, G j, add up to unity. 
Following the suggestion in ICappellaril d2002l) , we increased the 
minimal axial ratio of the Gaussians, qj, while not changing \ 2 
significantly, in order to make as large as possible the range of per- 
mitted inclinations. Parameters describing the parametrisation of 
the PSF and the deconvolved MGE model are presented in Tables|3] 
and|4] while a comparison of the combined ground and space-based 
light profiles and the MGE model is shown in Fig. [9] 

5.2 Simple dynamical models 

The basic assumption of our simple dynamical models is that 
clouds of ionised gas move in circular orbits within a thin disk. 
Th is is an often used app roach pioneered for black hole studies 
bv lMacchetto et all d 1997b and subseq uently applied to the i onised 
gas disk in Cen A by other authors (Ma rconi et al.ll200ll . 120061 ; 
Irlaring-Neumaver et alj200q) . 

The circular velocity is given by V C (R) = \J i?(d$/d_R), 
where the gravitational potential $ is as in eq. $2&. If mass-to-light 
ratio M/L is taken to be zero, the gravitational potential is described 
only by a 'point-mass' contribution ($.), and the circular velocity 




orcsec arcsec 



Figure 9. Comparison between the combined NICMOS/F222M and 
2MASS K-band photometry (open squares) with the MGE model (line) in 
seven angular sections as a function of radius. The right panels show the 
radial variation of the relative error along the profiles. 



Table 4. The parameters of the MGE model of the deconvolved K-band 
surface brightness of Cen A. 



j 


lj(L @ pc 2 ) 


Uj (arcsec) 




L J (xlO 9 L ) 


1 


70122.1720 


0.0506 


0.9500 


0.00030887 


2 


113639.0196 


0.4428 


1.0000 


0.040313 


3 


68662.2409 


1.2923 


1.0000 


0.20744 


i 


56428.6385 


2.9233 


1.0000 


0.87241 


5 


34809.6545 


5.3776 


1.0000 


1.8211 


6 


19386.9662 


13.6883 


0.9500 


6.2431 


7 


7184.3315 


37.4720 


1.0000 


18.250 


8 


3174.4729 


92.0318 


0.9823 


47.783 



simplifies to V C (R) = \J GM,/R. This model we call Kepler, 
while the more complex model, which includes the contribution of 
the stellar potential we call MGE. The calculation of the potential 
from the MGE parame trisation in that case is done as described in 
ICappellari etaild2002h . 

Thin disk models assume that the motion is purely circular 
and that the velocity dispersion of the gas is negligible. Our obser- 
vations of the Pa/3 disk are consistent with this assumption, but it 
breaks in the case of the [Fell] disk. We speculate that the increase 
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Table 5. The parameters of the best fit thin disk models 



model mass i 


PA 


M BH 


M/L 


x 2 


model [degr] 


[degr] 


[f0 7 M.] 


[K band] 




Pa/3 



pi 


MGE 


25 


-3 


q or; + 2.25 


0.0 + 1.3 


33.5 


p2 


MGE 


30 


-23 


45+ 1 - 


77.0 


p3 


MGE 


25 


-15 


9.0tl;g 


1^+1.55 


46.9 


pi 


MGE 


15 


-27 


3.5+1-0 


n 05 + O.85 


104.3 


p5 


Kepler 


25 


-3 


8.251?;?° 


0.0 


33.5 


p6 


Kepler 


30 


-23 


8.ol 2 2 ;° 


0.0 


78.1 


p7 


Kepler 


25 


-15 


9 5+ 20 


0.0 


47.0 


p8 


Kepler 


45 


-27 


5 o+ 13 

°- u -1.0 


0.0 


107.9 



Notes - We used two mass models: Kepler where the potential is given by 
the central point-mass, a mass model obtained from our MGE parametrisa- 
tion of the stellar light. There were 47 velocity data points used to constrain 
the models. 

in the velocity dispersion is related to the location of the [Fell] 
emission in the disk, the geometry of the warped disk and the tur- 
bulent medium with increasingly important [Fell] ionising shocks 
waves (see Section l531 for details). From this reasons, our observa- 
tions of [Fell] emission are not good tracers of the central potential 
and are not suitable for the determination of the black hole mass 
with the presented dynamical models. In what follows, we focus 
only on the Pa/3 emission, but we also construct a model for the 
[Fell] disk for comparison with observations (Figs[7jand[8}. 

To accurately reproduce the observed kinematics we include 
in the model the parametrisation of the emission-line surface 
brightness obtained by fitting a function to the measured flux, tak- 
ing the PA and the i of the gas disk into account. We tested differ- 
ent functional forms, finally choosing the following one which de- 
scribes the spatial distribution of the emission-lines surface bright- 
ness well: 



I = Io + h exp , 



(3) 



where ri is the scale radius and Io and Ii are the scale factors 
for the given gas component. The best fitting values for the Pa/3 
assuming for the orientation of the disk (PA — — 3°,i = 25°) 
are I = 1844.4, h = 32884.9 and n = 0'.'503. For other ori- 
entations we obtain slightly different values and use those in the 
corresponding dynamical models. 

When constructing dynamical models, we fix the geometry of 
the disk, use the above mentioned emission-line surface brightness 
parametrisation, and vary parameters defining the gravitational po- 
tential: M/L ratio and Al,. The resulting models of the circular 
velocity maps, convolve d with the PSF a nd the pixel size of the 
CIRAPSS observations dOian et alj[l995l) . are compared with the 
observed models. The best fitting model is determined by minimis- 
ing the x 2 — Ylj [(^wi>j — Vj)/AV,] 2 , similarly as done before. 



5.3 Results and discussion 

We calculated a range of dynamical models with different disk 
orientations, black hole masses and mass-to-light ratios. We ap- 
plied the disk orientation from Section |4~3l f or Pa/3 emission-line 
disk and also used values from the literature iMarcometalT 2006; 
lHaring-Neumaver et all2006t) . 

Table [5] presents the values of the best fitting models. Our 
best fitting orientation of the Pa/3 disk yields a model with M, = 




2.0X10 7 



Figure 10. A grid of models with different mass of the black hole (Mbh) 
and mass-to-light ratio (M/L) for the Pa/3 disk at i = 25° and PA = 
—3°. Over-plotted are the contours of the constant \ 2 . The first three con- 
tours correspond to 68.3%, 95.4% and 99.73% (thick contour) confidence 
regions for a distribution of 2 degrees of freedom. The best fitting model is 
shown with a diamond. 



(8.2514:25) x 10 7 M Q and M/L = 0.0±1.3. Uncertainties quoted 
for all presented models are at the 3<r level. The expected value 
of M, from the M, — a relation is excluded at a 3a confidence 
level, taking into accou nt the scatter of the relation of 0.34 dex 
dFerrarese & Fordl2005h . but our estimate, within its large uncer- 
tainties, is consisten t with other measurements in the literature at 
similar inclinations jMarconietaHl2006l :lH aring-Neumaver et al 



20061) . Values for M. obtained from stellar kinematics dSilge et al 
2005) are still higher (even their model for i — 20°) and formally 
excluded at a 3<r level. The estimated M/L ratio of the best fitting 
model is consistent with being zero confirming that within 1" the 
gravitational potential is almost entirely determined by the black 
hole and not dependant much on the parametrisation of the stellar 
light. 

Confidence levels over-plotted on the grid of different models 
with our best fitting disk orientation are presented in Fig. [Toj The 
reconstructed velocity map of the best model is shown in Fig. [7] 
and a cut along PA — 0° in Fig. [8] For completeness, we report 
here that using our best fitting orientation for the [Fell] disk we 
get Mbh = (O.llo °) X 10 7 A'7. and M/L = 2.0 ± 0.4. The 
untrustworthiness of this model can be seen in Figs. [7] (last panel) 
and[8](dash-dotted line). 

In the construction of the thin disk models the orientation of 
the disk is crucial, and the final parameters of the potential (M., 
M/L) will depend on it. We constructed other dynamical models 
(Table[5) using a set of disk geometries. For Pa/3: the best orienta- 
tion obtained from [Fell] kinemat ics (model p2), the best disk ori- 
entation from Ma rconi et al. d2006t) (model p3) and geometry of the 
best model fron^Haring-Neumav er et alj d2006h (model p4). Model 
p5 has the best Pa/3 disk geometry obtained from the kinematics, 
while models p6 - p8 have the same geometries as above (in the 
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Table 6. Additional disk models for given M. and M/L ratio. 



model 


mass 
model 


M B h 
[10 7 M.] 


M/L 
[K band] 


i 

[degr] 


PA 
[degr] 


x 2 


Pa/3 


P9 1 

pgi 

pgi 


MGE 
MGE 
MGE 
MGE 


6.1 
11.0 
24.0 
3.0 


0.7 
0.0 
0.68 
0.7 


25 ±5 
20 ±5 
20 ±1 
30±i° 


c + 10 
°-10 

oi? 3 

°-15 


37.3 
34.6 
128.6 
44.4 



same order), but these mass models include only the contribution 
from the black hole. 

In the last column of Table [5] x 2 values for each model are 
listed. The best model is pi, with a geometry determined from the 
kinematics derived in Section |4~31 Models pi and p5 are formally 
equally good, but other models (including p3 and p7) are excluded 
at the 3<7 level (for 2D distribution). The models are more sensitive 
to the inclination than to a change in PA. Models with Keplerian 
and MGE mass parametrisation give similar results, but the uncer- 
tainties in the estimated mass are somewhat smaller in Keplerian 
models, which have a fixed M/L. 

As an additional check of our results, we also constructed 
models with a given M, and M / L with the purpose to recover the 
measured disk orientations ( PA, i). We chose mo d els w ith param- 
eters from the literature: iHaring-Neumaver et al.l d2006h - model 
pgl. lMarconi et all d200r3) - model pg2. ISilge et alJd2005l) - model 
pg3. The models pg4 is a hybrid combining the value from the 
M. — a relation and the M/L ratio from other studies. Results 
are given in Table [6] 

As with previous models, most configurations give similar 
goodness of fit values, only model pg3 has a very different \ an d 
can be discarded. All models recover mutually consistent (within 
errors) geometrical parameters which are also in g ood agreement 
with v alues obtained in Section 14.31 Model pg2 ^Marconi et al.l 
d2006h geometry) is expectedly the closest t o our best fitting mode l 
(pi) in terms of \ 2 - Assuming M. from the lTremaine et al] J2002h 
relation gives a model that is close to be discarded with 3cr confi- 
dence level. 



6 CONCLUSIONS 

This paper presents observations of the active elliptical galaxy, 
NGC 5128 (Cen A), obtained with the integral-field spectrograph 
CIRPASS. In the spectra there are three emission-lines, [PII], [Fell] 
and Pa/3, with different spatial distributions and relative strengths. 
All lines are present in the central 1", [Fell] being the most ex- 
tended. Pa/? is also detected in two blobs, North and East of the 
nucleus. These blobs are in an excellent agreement with Paa blobs 
discovered by the HST imaging. Analysis and results are sum- 
marised in the following items. 

• The reconstructed total intensity image shows two regions 
with strong continuum contribution: in the centre and NE of the 
centre. The NE region overlaps with the position of the radio jet (in 
projection) and is located between the Pa/3 blobs. The region can 
also be associated with the knot Nl of the radio jet. This suggests a 
possible jet-induced star forming region in Cen A on scales of ~ 30 
pc, in addition to the well known large scale (lkpc) jet-induced star 
forming region. 



• The location of the East Pa/3 blob coincides with the kink in 
the radio jet suggesting a direct link between the two phenomena. 
If one assumes that the North Pa/3 blob was also induced by the jet, 
the position of the blobs show that the jet precessed for about 25°in 
projection. 

• We detect [PII] emission, which can be used together with 
[Fell] emission to discriminate between photo- and shock-induced 
ionisation of gas clouds. The [FeII]/[PII] ratio rises steadily with 
radius, being close, but above pure photo-ionisation levels. This 
suggests that shocks are more important with increasing distance 
from the centre. 

• We determined spatially resolved kinematics from Pa/3 and 
[Fell] emission-lines. Although the emission (of both lines) is di- 
vided in blobs or it is patchy at certain regions, it has a well defined 
and uniform sense of rotation and spatially separated emission re- 
gions are part of the same kinematic structure. The velocity dis- 
persion is low for Pa/3, but high for [Fell] emission with a plateau 
oriented at ~ —60° . The low dispersion is suggestive of a motion in 
a regular disk for Pa/3, but the plateau in the [Fell] velocity disper- 
sion is a clear evidence of a more complex distribution of ionised 
iron. 

• Kinematics of both emission-lines is remarkably similar in ap- 
pearance, as traced by ZVC: central regions are rotating around an 
EW axis, while at larger radii there is a twist in ZVC and the rota- 
tion is around NE - SW axis. The twist in ZVC suggest contribution 
from a bar-like perturbation or a warp in the disk on 20 pc scales. 

• We test the assumption that the emission in the central 1-2" 
is confined to a disk by constructing maps of circular velocity and 
comparing them with observations. This process is divided in two 
steps. We first determine the PA and (using this position angle) 
the inclination of the assumed disk. The best fitting parameters are: 
PA = -3° ± 10° and i = 25° ± 5 for Pa/3 and PA = -23° ± 3° 
and i = 30° ± 5 for [Fell]. Using these values we extract circu- 
lar velocity component and remove it from the data. The residual 
map of Pa/3 velocities is noisy without recognisable pattern and 
residuals are largely within measurement errors. [Fell] residuals are 
more complicated, showing a three-fold pattern and being some- 
what larger than measurement uncertainties. This can be traced to 
the twist in the ZVC. We conclude that the emission in the central 
1", the region which is likely to be within the sphere of influence of 
the central black hole, is consistent with being in circular motion in 
a disk, but outer regions show departures which could be attributed 
to perturbations (such as a bar) or to a warp extending inwards from 
1 kpc scales (as previously detected in Cen A) to 20 pc scales. 

• A shallower rotation curve of the [Fell] emission is likely re- 
lated to the high velocity dispersion and the geometry of the warped 
disk. We speculate that the excitation of [Fell] line happens also at 
larger radii from the central engine and the line-of-sight passes two 
times through the warped disk, such that at the spatial resolution of 
our observations the data contain a contributions of gas moving in 
two different directions (lowering the velocity, but increasing the 
velocity dispersion). This is supported by coinciding orientations 
of the high velocity dispersion region and the warped disk. Follow- 
ing this, we do not consider our [Fell] kinematics for determining 
the mass of the central black hole. 

• We construct simple dynamical models of the Pa/3 emission- 
line disk to determine the mass of the central black hole. For the 
geometry of disks we use kinematically determined orientation pa- 
rameters. We measure M. = 8.25i4;25 X 10 7 M . This value is 
consistent with previous measurements, while the predicted value 
of the black hole mass from the M. — a relation is excluded at the 
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3<t level (for a 2D distribution), taking into account the scatter of 
the relation. 

• We construct models with different geometrical orientations 
taken from previous studies, which give similar results, being 
marginally less good, but often statistically indistinguishable. We 
also construct models using literature values for M. and the value 
predicted by the M. — a relation, but leaving geometrical parame- 
ters free. These models confirm our original results. 
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